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ABSTRACT 

We describe the infrared properties of sources detected over ~ 36 deg 2 of sky in the GAMA 15-hr 
equatorial field, using data from both the Herschel Astrophysical Terahertz Large-Area Survey (H- 
ATLAS) and Wide-field Infrared Survey ( WISE). With 5<r point-source depths of 34 and 0.048 mJy 
at 250 /im and 3.4 /tin, respectively, we are able to identify 50.6% of the H- ATLAS sources in the 
WISE survey, corresponding to a surface density of ~ 630 deg -2 . Approximately two-thirds of these 
sources have measured spectroscopic or optical/near-IR photometric redshifts of z < 1. For sources 
with spectroscopic redshifts at z < 0.3, we find a linear correlation between the infrared luminosity 
at 3.4 fim and that at 250 pm, with ±50% scatter over ~ 1.5 orders of magnitude in luminosity, 

~10 9 — 10 10 - 5 L 0 . By contrast, the matched sources without previously measured redshifts (r > 20.5) 
have 250-350 pm flux density ratios that suggest either high-redshift galaxies (z > 1.5) or optically 
faint low-redshift galaxies with unusually low temperatures ( T < 20). Their small 3.4-250 pm flux 
ratios favor a high-redshift galaxy population, as only the most actively star-forming galaxies at low 
redshift (e.g., Arp 220) exhibit comparable flux density ratios. Furthermore, we find a relatively large 
AGN fraction (~30%) in a 12 pm flux-limited subsample of H-ATLAS sources, also consistent with 
there being a significant population of high-redshift sources in the no-redshift sample. 

Subject headings: surveys - cosmology: observations — galaxies: high-redshift - infrared radiation - 
galaxies: statistics - galaxies: general 
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1. INTRODUCTION 

The advent of submillimeter astronomy has opened a 
new window into the Universe, allowing us to probe dusty 
star-forming galaxies at high redshift, as well as cold 
dust in nearby galaxies. Along with the recent success 
of the balloon-borne BLAST experiment (Pascale et al. 
2008; Eales et al. 2009), the launch of the Herschel Space 
Observatory 1 (Pilbratt et al. 2010) allows us to probe 
wavelengths from 55 to 672 pm to a spatial resolution of 
< 10". Its largest open-time key project, the Herschel 
Astrophysical Terahertz Large-Area Survey (H-ATLAS, 
Eales et al. 2010a), will observe ~ 550 deg 1 2 of sky and 
detect more than 300,000 galaxies. 

Based upon a preliminary cross-identification with the 
Galaxy and Mass Assembly survey (GAMA, Driver et al. 
2011) and Sloan Digital Sky Survey (SDSS, Abazajian 
et al. 2009), ~50% of the sources detected in the H- 
ATLAS survey are at z < 1 (Smith et al. 2011). In 
addition, an analysis of the FIR colors of the remaining 
H-ATLAS sources with flux densities above 35 mJy at 
350 pm and > 3tr detections at 250 pm and 500 pm sug- 
gests an average redshift of z ~ 2 (Amblard et al. 2010; 
Lapi et al. 2011). 

Another method of constraining the redshift distribu- 
tion of FIR-selected sources is to look for counterparts 
in the near and mid infrared (NIR and MIR). Fleuren 
et al. (in prep) have performed source matching to sur- 
vey data from the VISTA Kilo-degree Infrared Galaxy 
survey (VIKING, Sutherland et al., in prep), but these 
data only extend to ~ 2 pm. Starting in December 2009, 
the Wide-field Infrared Survey Explorer 2 ( WLSE, Wright 
et al. 2010) began its mission to observe the entire sky 
in four bands, ranging from 3.4 to 22 pm, at < 12" res- 
olution. Following the first public data release in April 
2011, ~ 24,000 deg 2 of IR images and source catalogs be- 
came available to the public, including 36 deg 2 of sky in 
the equatorial plane covered by the H-ATLAS survey. 

The purpose of this letter is to describe the infrared 
properties of identified WISE counterparts to H-ATLAS 
sources within the GAMA 15-hr field (G15). Through- 
out we will assume a concordance cosmology with Ho = 
71 km s -1 Mpc -1 , f2 m = 0.27, and Oa = 0.73 (Spergel 
et al. 2007). In a subsequent paper, we will report WISE 
cross-identifications over the entire H-ATLAS area and 
perform spectral energy distribution (SED) fits to the 
matched sources. 

2. DATA AND METHODOLOGY 

Herschel observations of the G15 field include imaging 
data at 250, 350, and 500 pm from the SPIRE instru- 
ment. with respective beam FWHM of 18"1, 25"2, and 
36" 6 (Griffin et al. 2010), as well as 100 and 160pm imag- 
ing from the PACS instrument (Poglitsch et al. 2010). 
The field subtends approximately 12° in right ascension 
and 3° in declination and is centered on the equatorial 
plane at a = 14 h 30 m. The H-ATLAS G15 source cat- 
alogue (Dunne et al., in prep) is constructed using the 
same method as the Science Demonstration Phase cat- 

1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with 
important participation from NASA. 

2 For a description of the WISE mission, see 

http://wise2.ipac.caltech.edu/docs/release/prelim/expsup/ 


alogue (Rigby et al. 2011). It contains 27,481 sources 
detected at > 5<r in any of the three SPIRE bands and 
reaches point-source depths of 34, 40, and 44 mJy at 250, 
350, and 500 pm, respectively. 

The WISE first public data release contains four-band 
coverage of the entire G15 field to 5<r point-source depths 
of 0.048, 0.10, 0.73, and 5.9 mJy at 3.4, 4.6, 12, and 
22 pm. The angular resolution in these bands is 6"1, 6"4, 
6"5, and 12"0, respectively (Wright et al. 2010). There 
are ~ 240,000 IVJS'E sources within G15 detected at > Icr 
in at least one of the four WISE bands, corresponding to 
~ 0.14 WISE sources per 10"-radius aperture. In the 
WISE preliminary release catalog, close pairs of sources 
are not deblended for separations <9", so when a back- 
ground or foreground source appears near an H-ATLAS 
source position, it will often be blended with the true 
counterpart. 

In addition to the infrared data from WISE and H- 
ATLAS, the G15 region has spectroscopic redshifts from 
the GAMA survey (r < 19.4, Driver et al. 2011), and 
photometric redshifts obtained using optical/NIR pho- 
tometry from SDSS, VIKING, and the UK Infrared Deep 
Sky Survey Large Area Survey (UKIDSS-LAS, Lawrence 
et al. 2007). Of the 27,481 H-ATLAS sources in the re- 
gion, 16% have reliable spectroscopic redshifts and 39% 
have photometric redshifts. 

As a result of the high space density of 3.4 pm sources, 
a naive matching to the H-ATLAS source positions will 
result in a non-negligible number of misidentifications. 
A significant fraction of these misidentifications will be 
foreground stars that are blue in all WISE bands and 
undetectable in H-ATLAS, so we can reduce the con- 
tamination rate of our matched catalog by considering 
only WISE sources with [3.4] — [4.6] > 0. The magnitude 
distribution of all background/foreground WISE sources 
is such that ~ 20% are bluer than this limit, as compared 
to < 1% of sources within 10" of an H-ATLAS position. 

To estimate the fraction of H-ATLAS sources with 
detectable counterparts in the remainder of the WISE 
source catalog, we use the method of Fleuren et al., in 
prep, where it is shown that by counting non-matches 
within a fixed distance from the H-ATLAS positions, we 
can estimate the detection rate, 

Qo = 1 - ( 1 ) 

where S is the fraction of unmatched positions in the 
H-ATLAS catalog and B is the fraction of unmatched 
random positions. This expression holds under the con- 
dition that there are a negligible number of true WISE 
matches outside of the matching radius. Matching all 
WISE sources within 10" (within which we expect > 99% 
of the true matches to lie) and using 10 5 random posi- 
tions within the G15 field, we find Qo = 0.632 ± 0.004. 
This is larger than the value found in the SDSS galaxy- 
catalog (Qo = 0.583, Smith et al. 2011), but smaller than 
in the VIKING Ks band (Qo = 0.75, Fleuren et al., in 
prep). By contrast, we find Qo = 0.012±0.002 within the 
subset of WISE sources with [3.4] — [4.6] < 0. suggesting 
that this color cut was effective in removing stars. 

In Fig. 1, we use Eq. 1 to estimate the IR detection rate 
as a function of FIR flux density. We find that WISE 
3.4 pm detection rates are > 95% for sources brighter 
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Fig. 1. — NIR and MIR detection rates plotted as a function 
of limiting FIR flux density at 250pm (solid), 350pm (dashed), 
and 500pm (dotted). Detection rates are shown for 3.4pm (upper 
left panel), 4.6pm (upper right), 12 pm (lower left), and 22 pm 
(lower right) sources, where the corresponding 5<r detection limits 
are 0.06, 0.9, and 3.6 mJy, respectively. Note that the large sky 
density of 3.4 pm sources can lead to ambiguity in the identification 
of some sources, so identification rates are typically ~ 10-20% lower 
than the corresponding detection rates. 


than 100 mJy at 250 pm. The vast majority of these 
sources (78%) have measured spectroscopic or photomet- 
ric redshifts < 0.5, suggesting that the low-redshift mode 
of star-forming galaxies dominates sub-mm sources above 
this flux density. Below this FIE flux density, the NIR 
and MIR detection rates drop rapidly as the high-redshift 
mode accounts for an increasing fraction of the sub-mm 
sources (see Fig. 2). 

By contrast, ~ 50% of sources brighter than 100 mJy 
at 500 pm are expected to be lensed z > 1 galaxies (Ne- 
grello et al. 2010; Hopwood et al. 2011). Of the 46 
such sources in the G15 region, — 20% are undetected 
at 12 and 22 pm, most likely because the 500 pm band- 
pass shifts higher in the blackbody curve at high redshift, 
while the rest-frame MIR flux density declines blueward 
of 22 pm (e.g., Rieke et al. 2009). 

The identification of individual NIR and MIR counter- 
parts is more subtle, as we want to eliminate as many of 
the false matches as possible. Here, we use the likelihood 
ratio technique of Sutherland k Saunders (1992), which 
was implemented for the H- ATLAS survey in Smith et al. 
(2011, hereafter Sll). Because > 99.9% of sources in the 
WISE catalog are detected at 3.4 pm at > 7<r, we use 
this band for source identification. The likelihood that 
a given WISE source is a counterpart to an H-ATLAS 
source is a function of the radial probability distribu- 
tion, f(r), the 3.4pm magnitude distribution of non- 
counterparts, n(m), and the 3.4 pm magnitude distribu- 
tion of detected H-ATLAS counterparts, q(m): 


Q(m)f(r) 

n(m) 


where f(r) is given by 


Hr) = 
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-exp 


r 2 \ 

hz ) 


( 2 ) 


Fig. 2. — Distribution of the 3.4— 250/iim flux density ratio as a 
function of 250 fun flux density for H-ATLAS sources with WISE 
3.4 fun IDs. Individual points are plotted for sources with spec- 
troscopic redshifts from the GAMA survey, with the color of the 
point indicating the redshift. The sources without either photo- 
metric or spectroscopic redshifts are indicated by black points. 
Dashed lines indicate the approximate limits of the H-ATLAS and 
WISE surveys. The majority of sources with spectroscopic red- 
shifts (r < 19.4) are well detected by WISE, while the no-redshift 
sources have small flux density ratios that are inconsistent with the 
majority of “normal” low-redshift galaxies. 


We measure n(m) directly from the WISE prelim- 
inary release catalog, while we determine q{m) by 
first measuring the magnitude distribution of 3.4 pm 
sources within 10" of the H-ATLAS source positions 
and then subtracting the magnitude distribution of 
background/ foreground sources normalized by the area 
within this radius. 

Following Sll, we assume the astrometric uncertainty 
to be circularly symmetric and to have a dependence on 
the 250 pm signal-to-noise ratio (SNR) given by, 


Up os 


0.655 


FWHM 
SNR ’ 


( 4 ) 


where FWHM= 18.1" is the Hers chel beam width. While 
the WISE random astrometric uncertainties are negli- 
gible (~0'. / 3), the WISE positions in the preliminary 
release catalog are known to suffer from highly non- 
Gaussian systematic offsets that can be as large as 1", so 
we apply a lower limit of cr pos > 1.5". 

Once we have computed the likelihood, Lj, for a pos- 
sible counterpart j, we determine the reliability by sum- 
ming over all possible counterparts, 




Lj 

E» £t + (i-Qo)‘ 


( 5 ) 


We consider any 3.4 pm source with R > 0.8 to be a 
reliable counterpart. After performing this procedure 
separately for WISE sources with star- and galaxy-like 
[3.4] — [4.6] colors, we obtain a combined total of 13,898 
WISE counterparts, or 50.6 ±0.4% of the H-ATLAS cat- 
alog. The expected number of false matches can be ob- 
tained with 

A'falselD — y~/l — lit)- 
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( 3 ) 
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Fig. 3. — FIR flux density ratio distribution for 311 H-ATLAS 
sources with WISE identifications, but without measured redshifts. 
Colored points indicate the median redshift of a set of modified 
blackbody models with temperatures uniformly distributed be- 
tween 10 and 60 K and emissivity parameters between 0 and 2. 
Dashed lines indicate the median flux density ratios of the sources 
lacking redshifts. The majority of the sources without redshifts ap- 
pear to be at z ~ 2. Inset. Combined distribution of spectroscopic 
and photometric redshifts for H-ATLAS sources with detections at 
3.4/rm (solid histogram), 12 fim (dashed), and 22 /rm (hatched). 
There are 73 sources at z > 1 not shown. 

We estimate that there are 369 false identifications in 
our sample, corresponding to a contamination rate of 
2.7%. The majority of our IDs have galaxy-like colors, 
but we were able to identify 14 objects with star-like 
colors ([3.4] — [4.6] < 0) that weren’t already identified 
as galaxies with SDSS matching. This corresponds to an 
upper limit of 0.05% for the fraction of H-ATLAS sources 
identified as stars and detected in WISE. 

Comparing the individual WISE identifications to 
SDSS cross-identifications, we find that 84.7% of the 
10, 709 G15 H-ATLAS sources with R > 0.8 SDSS identi- 
fications also had R > 0.8 WISE identifications. The me- 
dian separation between the WISE and SDSS sources is 
0 , /6, consistent with the astrometric uncertainties in the 
current version of the WISE catalog. In addition, ~ 6% 
of these sourcess have separations > 3", consistent with 
the expected contamination rates of the WISE (~ 3%) 
and SDSS (—5%) matched samples. 

3. DEMOGRAPHICS 

The redshift distribution of 250 pm-selected sources 
is likely bimodal, with populations of both moderately- 
star-forming galaxies at 2 < 1 and high-redshift starburst 
galaxies at 2 ~ 2 (Smith et al. 2011; Lapi et al. 2011). 
Because the latter are typically faint in the observed- 
frame optical and UV, the vast majority of objects with 
measured spectroscopic or photometric redshifts are in 
the low-redshift. mode. Of the H-ATLAS sources with 
spectroscopic and photometric redshifts (see the inset of 
Fig. 3), we successfully identify 85.6% in WISE. 

Although most of the WISE identifications with known 
redshifts are at 2 < 0.8 (> 98%), such a sample is biased 
toward low redshift because redshift measurements re- 
quire detectable flux in the UV, optical, or near-infrared. 
For the remaining sources, we can estimate redshifts us- 


ing their FIR flux density ratios (e.g., Amblard et al. 
2010; Schulz et al. 2010), which give an indication of the 
location of the peak of the FIR dust emission. In Fig. 3, 
we compare the SPIRE colors of our matched catalog to 
a suite of 10 6 modified blackbody spectral energy distri- 
butions (including uncertainties) with temperatures uni- 
formly distributed between 10 and 60 K and emissivitiy 
parameters between 0 and 2 (Amblard et al. 2010). We 
indicate on the plot the median flux density ratios of the 
4, 677 H-ATLAS sources with WISE identifications but 
no redshifts ( “no-redshift” sample, dashed lines), as well 
as the individual flux density ratios of 311 such sources 
with > 3a detections in all three SPIRE bands. 

The no-redshift sample has median FIR colors suggest- 
ing 2 ~ 2, while the subset with three-band SPIRE de- 
tections may have even higher redshifts, perhaps as large 
as 2 ~ 3.5. The redshifts derived in this way are highly 
uncertain, as redshift is degenerate with dust temper- 
ature in the FIR and it is possible that these galaxies 
have typical dust temperatures different than 35 K. How- 
ever, in order for the no-redshift sample to be at 2 ~ 0, 
they would need to have very cold dust temperatures 
(T ~ 12 K) in addition to being very faint in the optical 
(r > 20.5). Although temperatures as low as 10 K have 
been seen in H-ATLAS galaxies at 2 < 1, T ~ 25 — 30 K 
is more typical (Dye et al. 2010). 

Further evidence that the no-redshift sources are at 
z> 0 can be found in the 3.4-250 pm ratios (see Fig. 2). 
H-ATLAS sources with GAMA redshifts have NIR-FIR 
flux density ratios that decrease from ~ 0.009 at 2 ~ 0 to 
~ 0.006 at z ~ 0.35. This trend, which continues toward 
higher redshift, is due primarily to the larger fe-correction 
at 3.5 /am than at 250 pm (e.g., Rieke et al. 2009). The 
median flux density ratio for the no-redshift sources, by 
comparison, is 0.0023. Although low-redshift galaxies 
will occasionally exhibit such small NIR-FIR flux density 
ratios - Arp 220, for example, has / 3 . 4//250 = 0.0019 - 
such objects are rare, actively star-forming, and unlikely 
to have dust temperatures as small as ~ 15 K. 

We can shed even further light on the properties of 
the no-redshift sample by examining their position in 
WISE color space (see Fig. 4). The technique is simi- 
lar to those developed for a set of IRAC filters at com- 
parable wavelengths (Sajina et al. 2005; Stern et al. 
2005). It was demonstrated that the diagram can be 
divided into stars and early-type galaxies (lower left cor- 
ner, [4.6] — [12] < 1), star-forming galaxies (lower right 
corner), and active galactic nuclei (AGNs), where AGNs 
are selected largely based on their power-law emission in 
the NIR and MIR. This was adapted for WISE by Jarrett 
et. al. (2011), 

[4.6] - [12] > 2.2 

[4.6] - [12] < 4.2 

[3.4] - [4.6] < 1.7 (7) 

[3.4] - [4.6] > 0.1 ([4.6] - [12]) + 0.38. 

Of the 4,959 H-ATLAS sources with WISE identifica- 
tions, 2 < 0.5, and detections in 3.4, 4.6, and 12 pm, the 
majority are star-forming galaxies, with < 0.2% main 
sequence stars (which is consistent with the SDSS/H- 
ATLAS matching done by Thompson et al. 2010), and a 
small AGN fraction (0.057), 
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Fig. 4. — IR color distribution of ATLAS- WISE matches with 
three-band WISE photometry, including sources with z < 0.5 
(color coded by redshift) and no-redshift sources (black points). 
The diagram can be divided into stars (lower left corner), star- 
forming galaxies (lower right corner), and AGNs (dashed region). 
The color distribution of all three-band WISE detections (includ- 
ing sources unmatched to H- ATLAS) is indicated by contours. The 
majority of identified sources are normal star-forming galaxies, 
with a small fraction of AGNs. Sources without measured red- 
shifts tend to be redder in [3.4] — [4.6] and have a larger AGN 
fraction than sources with spectroscopic or photometric redshift 
measurements, z < 0.5. 

For the 9% of no-redshift sources with three-band 
WISE photometry, we find a much larger AGN fraction 
(0.30), but are subject to a selection bias toward AGNs 
due to their additional flux from warm dust at 12 pm 
(Hainline et al. 2010). Coppin et al. (2010) performed 
SED fits on a set of blank-field SMG samples and esti- 
mated that ~ 15% of SMGs are dominated by an AGN 
in the mid-infrared. 

4. FIR-NIR PROPERTIES 

Emission from dust in the FIR is often used as an in- 
dicator of the total star formation rate in galaxies (e.g., 
Kennicutt 1998), under the assumption that most of the 
dust heating is provided by young stars and that the star- 
forming regions are optically thick. The 250 gm band- 
pass samples the cooler dust emission, which can arise 
from both star-forming clouds and diffuse regions in the 
ISM (e.g., Eales et al. 2010b; Dunne et al. 2011). By 
contrast, 3.4 pm emission is dominated by stellar sources 
and is a probe of the stellar mass for galaxies without a 
strong AGN contribution (e.g., Stern et al. 2005). 

In Fig. 5, we plot 3.4 pm luminosity as a function of 
250 pm luminosity for H- ATLAS sources with spectro- 
scopic redshifts, 0.05 < 2 < 0.3. We fe-correct 3.4 pm 
using a power law interpolation of the [3.4] — [4.6] color. 
For the 250 pm flux densities, the majority of the galax- 
ies in our sample do not have sufficient FIR flux density 
to obtain reliable dust temperatures, as is needed for a 
proper k correction. Instead, we use the median dust 
temperature of 26 K found by Dye et al. (2010) using 
PACS and SPIRE data from the H- ATLAS SDP of 1,346 
0.1 < z < 0.5 sources. Although the r < 19.4 limit of 
the GAMA spectroscopic survey restricts the range of 
luminosities that can appear on the plot, we estimate 



Fig. 5. — Luminosity at 3.4 pm as a function of luminosity at 
250 p m for H- ATLAS sources with spectroscopic redshifts, 0.05 < 
z < 0.3. Both luminosities are rest-frame fe-corrected. Small 
squares are colored according to their redshift. Red stars indi- 
cate AGN, selected with the criteria given in Eq. 7. Dashed lines 
indicate the approximate selection limits for sources at z = 0.05 
and 0.2. The majority of the AGN lie above the relationship, sug- 
gesting that 1 , 3.4 contains a non-negligible contribution from hot 
dust emitting near the central black hole. 

that GAMA misses only ~ 5% of matched 2 < 0.3 galax- 
ies based upon the deeper (r < 20.5) subsample of H- 
ATLAS galaxies with photometric redshifts. 

We find that the majority of the IR-color-selected AGN 
(red points) lie above the relationship for star-forming 
galaxies, as T 3.4 will contain a non-negligible contribu- 
tion from hot dust emitting near the central black hole 
in addition to emission from low-mass stars. However, if 
we exclude AGNs, we find an approximately linear cor- 
relation, with a best-fit power law index, a = 0.98lo(o5- 
The quoted systematic uncertainties on the power law 
index were determined by allowing for a range of possi- 
ble dust temperatures, 18 < T < 34 K, when performing 
250 pm ^-corrections. The intrinsic scatter about this 
relationship is 0.18 ± 0.01 dex, or ~ 50%. 

The existence of a linear correlation between L 2 50 and 
L3.4 across one and a half decades in luminosity suggests 
a close relationship between the cold dust probed by the 
FIR and the stellar mass probed by the NIR. Previous 
indications with Herschel have shown the 250 pm lumi- 
nosity density to tightly correlate with both the 24 pm 
luminosity density, a star formation rate indicator, and 
the total infrared luminosity (Elbaz et al. 2010). This 
fact, coupled with a linear correlation between stellar 
mass and star formation rate for the general star-forming 
galaxy population (Daddi et al. 2007; Elbaz et al. 2007; 
Noeske et al. 2007, Donoso et al., in prep), suggests that 
the majority of the low-redshift H-ATLAS sources are 
actually “normal” star-forming galaxies. 

This publication makes use of data products from the 
Wide-field Infrared Survey Explorer, which is a joint 
project of the University of California, Los Angeles, 
and the Jet Propulsion Laboratory/ California Institute 
of Technology, funded by the National Aeronautics and 
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Space Administration. 

Herschel is an ESA space observatory with science in- 
struments provided by European-led Principal Investiga- 
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